Search for: All records

The engulfment of substellar bodies (SBs), such as brown dwarfs and planets, by giant stars is a possible explanation for rapidly rotating giants, lithium-rich giants, and the presence of SBs in close orbits around subdwarfs and white dwarfs. We perform three-dimensional hydrodynamical simulations of the flow in the vicinity of an engulfed SB. We model the SB as a rigid body with a reflective surface because it cannot accrete. This reflective boundary changes the flow morphology to resemble that of engulfed compact objects with outflows. We measure the drag coefficients for the ram-pressure and gravitational drag forces acting on the SB, and use them to integrate its trajectory inside the star. We find that engulfment can increase the luminosity of a 1M_⊙star by up to a few orders of magnitude. The time for the star to return to its original luminosity is up to a few thousand years when the star has evolved to ≈10R_⊙and up to a few decades at the tip of the red giant branch (RGB). No SBs can eject the envelope of a 1M_⊙star before it evolves to ≈10R_⊙if the orbit of the SB is the only energy source contributing to the ejection. In contrast, SBs as small as ≈10M_Jupcan eject the envelope at the tip of the RGB. The numerical framework we introduce here can be used to study planetary engulfment in a simplified setting that captures the physics of the flow at the scale of the SB.

 

Common-envelope evolution is a stage in binary system evolution in which a giant star engulfs a companion. The standard energy formalism is an analytical framework to estimate the amount of energy transferred from the companion's shrinking orbit into the envelope of the star that engulfed it. We show analytically that this energy transfer is larger than predicted by the standard formalism. As the orbit of the companion shrinks, the mass it encloses becomes smaller, and the companion is less bound than if the enclosed mass had remained constant. Therefore, more energy must be transferred to the envelope for the orbit to shrink further. We derive a revised energy formalism that accounts for this effect, and discuss its consequences in two contexts: the formation of neutron star binaries, and the engulfment of planets and brown dwarfs by their host stars. The companion mass required to eject the stellar envelope is smaller by up to 50% , leading to differences in common-envelope evolution outcomes. The energy deposition in the outer envelope of the star, which is related to the transient luminosity and duration, is up to a factor of ≈7 higher. Common-envelope efficiency values above unity, as defined in the literature, are thus not necessarily unphysical, and result at least partly from an incomplete description of the energy deposition. The revised energy formalism presented here can improve our understanding of stellar merger and common-envelope observations and simulations. 

The engulfment of substellar bodies (SBs) such as brown dwarfs and planets has been invoked as a possible explanation for the presence of SBs orbiting subdwarfs and white dwarfs, rapidly rotating giants, and lithium-rich giants. We perform three-dimensional hydrodynamical simulations of the flow in the vicinity of an SB engulfed in a stellar envelope. We model the SB as a rigid body with a reflective boundary because it cannot accrete. This reflective boundary changes the flow morphology to resemble that of engulfed compact objects with outflows. We measure the drag coefficients for the ram pressure and gravitational drag forces acting on the SB, and use them to integrate its trajectory during engulfment. We find that SB engulfment can increase the stellar luminosity of a 1M⊙ star by up to a few orders of magnitude for timescales of up to a few thousand years when the star is ≈10R⊙ and up to a few decades at the tip of the red giant branch. We find that no SBs can eject the envelope of a 1M⊙ star before it evolves to ≈10R⊙ . In contrast, SBs as small as ≈10MJup can eject the envelope at the tip of the red giant branch, shrinking their orbits by several orders of magnitude in the process. The numerical framework we introduce here can be used to study the dynamics of planetary engulfment in a simplified setting that captures the physics of the flow at the scale of the SB. 

The engulfment of substellar bodies (SBs) such as brown dwarfs and planets has been invoked as a possible explanation for the presence of SBs orbiting subdwarfs and white dwarfs, rapidly rotating giants, and lithium-rich giants. We perform three-dimensional hydrodynamical simulations of the flow in the vicinity of an SB engulfed in a stellar envelope. We model the SB as a rigid body with a reflective boundary because it cannot accrete. This reflective boundary changes the flow morphology to resemble that of engulfed compact objects with outflows. We measure the drag coefficients for the ram pressure and gravitational drag forces acting on the SB, and use them to integrate its trajectory during engulfment. We find that SB engulfment can increase the stellar luminosity of a 1M⊙ star by up to a few orders of magnitude for timescales of up to a few thousand years when the star is ≈10R⊙ and up to a few decades at the tip of the red giant branch. We find that no SBs can eject the envelope of a 1M⊙ star before it evolves to ≈10R⊙ . In contrast, SBs as small as ≈10MJup can eject the envelope at the tip of the red giant branch, shrinking their orbits by several orders of magnitude in the process. The numerical framework we introduce here can be used to study the dynamics of planetary engulfment in a simplified setting that captures the physics of the flow at the scale of the SB. 

The coalescence of two neutron stars was recently observed in a multi-messenger detection of gravitational wave (GW) and electromagnetic (EM) radiation. Binary neutron stars that merge within a Hubble time, as well as many other compact binaries, are expected to form via common envelope evolution. Yet five decades of research on common envelope evolution have not yet resulted in a satisfactory understanding of the multi-spatial multi-timescale evolution for the systems that lead to compact binaries. In this paper, we report on the first successful simulations of common envelope ejection leading to binary neutron star formation in 3D hydrodynamics. We simulate the dynamical inspiral phase of the interaction between a 12 M⊙ red supergiant and a 1.4 M⊙ neutron star for different initial separations and initial conditions. For all of our simulations, we find complete envelope ejection and a final orbital separation of ≈1.1 - 2.8R⊙ , leading to a binary neutron star that will merge within 0.01-1 Gyr. We find an αCE -equivalent efficiency of ≈0.1 - 0.4 for the models we study, but this may be specific for these extended progenitors. We fully resolve the core of the star to ≲0.005R⊙ and our 3D hydrodynamics simulations are informed by an adjusted 1D analytic energy formalism and a 2D kinematics study in order to overcome the prohibitive computational cost of simulating these systems. The framework we develop in this paper can be used to simulate a wide variety of interactions between stars, from stellar mergers to common envelope episodes leading to GW sources.